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GROUND STATES OF LARGE BOSONIC SYSTEMS: THE 
GROSS-PITAEVSKII LIMIT REVISITED 

PHAN THAnH NAM, NICOLAS ROUGERIE, AND ROBERT SEIRINGER 


Abstract. We study the ground state of a dilute Bose gas in a scaling 
limit where the Gross-Pitaevskii functional emerges. This is a repulsive 
non-linear Schrodinger functional whose quartic term is proportional 
to the scattering length of the interparticle interaction potential. We 
propose a new derivation of this limit problem, with a method that 
bypasses some of the technical difficulties that previous derivations had 
to face. The new method is based on a combination of Dyson’s lemma, 
the quantum de Finetti theorem and a second moment estimate for 
ground states of the effective Dyson Hamiltonian. It applies equally 
well to the case where magnetic fields or rotation are present. 
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1. Introduction 

The rigorous derivation of effective non-linear theories from many-body 
quantum mechanics has been studied extensively in recent years, motivated 
in part by experiments in cold atom physics. For bosons, the emergence 
of the limit theories can be interpreted as due to most of the particles oc¬ 
cupying the same quantum state: this is the Bose-Einstein condensation 
phenomenon, observed first in dilute alkali vapors some twenty years ago. 

The parameter regime most relevant for the description of the actual phys¬ 
ical setup is the Gross-Pitaevskii limit. It is also the most mathematically 
demanding regime considered in the literature so far, see [znuziiiH] for the 
derivation of equilibrium states and [iniEllIlEl! for dynamics (more exten¬ 
sive lists of references may be found in [201 [25l [2]). The main reason for 
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this sophistication is the fact that interparticle correlations due to two-body 
scattering play a leading order role in this regime. The goal of this paper 
is to present a method for the derivation of Gross-Pitaevskii theory at the 
level of the ground state that is conceptually and technically simpler than 
existing proofs, in particular that of [18] which was so far the only method 
applicable when an external magnetic field is present. 

Our setting is as follows: we consider N interacting bosons in the three- 
dimensional space described by the many-body Schrodinger Hamiltonian 

N 

HN = '^hj+ ^ WNiXj-Xk) (1.1) 

j=l l^j<k^N 

acting on the space permutation-symmetric square 

integrable functions. The one-body operator is given by 

h:= (-iV + H(x))2 + H(x) 

with a magnetic (or a rotation) field A satisfying 

H e Lf„^(M^M^), lim |H(a:)|e-'’l^l = 0 (1.2) 

|ir|—>-oo 

for some constant 6 > 0 and an external potential V satisfying 

0 ^ H G Lj'Qj,(M^), lim H(x) = -|-oo. (1.3) 

|x|—>-oo 

We thus consider non-relativistic particles in a trapping potential, possibly 
under the influence of an effective magnetic field, which might be due to 
rotation of the sample or the interaction with optical fields. 

The particles interact pairwise via a repulsive potential wn given by 

WNix) = N^w{Nx), (1.4) 


where w is a fixed function which is non-negative, radial and of finite range, 
i.e., l(lx| > Ro)w{x) = 0 for some constant Rq > 0. Different scalings of 
the interaction potential have been considered in the literature, of the form 

W/3,n{x) = ^N^^w{N>^x) (1.5) 

with 0^/3^ 1. The prefactor makes the interaction energy in (|l.lh of 
the same order as the one-particle energy. Indeed, if /3 > 0, then 

N^^w{N^x) iv^ 6o ( 1 . 6 ) 

weakly and thus the interaction potential should be thought of as 

leading to a bounded interaction energy per pair of particles. Generally 
speaking, the larger the parameter /3, the faster the potential converges to 
a point interaction, and thus the harder the analysis. Note that the cases 
/3 < 1/3 and /3 > 1/3 correspond to two physically rather different scenarios; 
in the former the range of the potential is much larger than the typical inter¬ 
particle distance and we should expect many weak collisions; while in 

the latter we rather have very few but very strong collisions. In this paper 
we consider the most interesting case j3 = 1 where the naive approximation 
(|1.6I) does not capture the leading order behavior of the physical system. 
In fact, the strong correlations at short distances 0{N~^) yield a nonlinear 
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correction, which essentially amounts to replacing the coupling constant f w 
by (87r)x the scattering length of w. 

Let us quickly recall the definition of the scattering length; a more com¬ 
plete discussion can be found in [20l Appendix C]. Under our assumption 
on w, the zero-energy scattering equation 

{—2A + w{x))f{x), lim /(x) = 1, 

|x|—>-oo 

has a unique solution and it satisfies 

fix) = 1 - > ^0 

\x\ 

for some constant a ^ 0 which is called the scattering length of w. In 
particular, if w is the potential for hard spheres, namely w{x) = -|-oo when 
|x| < Rq and rc(x) = 0 when |x| ^ then the scattering length of w 
is exactly Rq. In a dilute gas, the scattering length can be interpreted 
as an effective range of the interaction: a quantum particle far from the 
others is felt by them as a hard sphere of radius a. A useful variational 
characterization of a is as follows: 

Sira = inf i [ 2|V/p-h u;|/p, lim /(x) = . (1.7) 

UR3 |x|^-oo J 

Consequently, Svra is smaller than f w (the strict inequality can be seen by 
taking the trial function 1 — Xg with g G C^(M^, M) satisfying g{x) = 1 when 
|x| < Rq, and A > 0 sufficiently small). Moreover, a simple scaling shows 
that the scattering length of wjy = N‘^w{N.) is a/N. 

We are going to prove that the ground state energy and ground states of 
Hj\f converge to those of the Gross-Pitaevskii functional 

Scpiu) := {u,hu) -|-dvra / \u{x)\^dx (1-8) 

in a suitable sense. Note that the occurrence of the scattering length in 
(jl.Sp is subtle: this functional is not obtained by testing Hn with factorized 
states of the form (which would lead to a functional with dvro replaced 
by (1/2) f w). Taking into account the short-range correlation structure 
which gives rise to (II.8h is the main difficulty in the proof of the following 
theorem, which is our main result. 


Theorem 1.1 (Derivation of the Gross-Pitaevskii functional). 

Under conditions (IL21), (1131) and dni), we have 


lim inf 


(T,77v^) 

N 


inf TGp('ii) =: ecp. 
l“lh2(R3) = l 


(1.9) 


Moreover, if n is an approximate ground state for Hp 4 , namely 


(Tv, LIv'I'Ar) 

hm -—- 

V—>-CXD N 


fiGP, 


then there exists a subsequence T and a Borel probability measure // sup¬ 
ported on the set of minimizers of £gp{u) such that 


lim Tr 
1^00 

- j 

f \u^’^){u^’^\d^i{u) 

= 0, VA: G N 


( 1 . 10 ) 
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where 7 ^^^ = Trfc_|_i_j. 7 v |'I^Ar)(^Af| is the k-partiele reduced density matrix of 
'I'Ar. In particular, if£Qp{u) subject to \\u\\i 2 = 1 has a unique minimizeruo 
has a unique minimizer uq (up to a eomplex phase), then there is complete 
Bose-Einstein condensation 


lim Tr 

N^OO 




\/k G N. 


( 1 . 11 ) 


The energy upper bound in (II.bp was proved in [211 [25] (see also [a Ap¬ 
pendix A] for an alternative approach). The energy lower bound in (|1.9p and 
the convergence of one-particle density matrices were proved in [18]. The 
simpler case A = 0 had been treated before in [21] (ground state energy) and 
[17] (condensation). In this case, the uniqueness of the Gross-Pitaevskii min¬ 
imizer Uq follow from a simple convexity argument. The result in Theorem 
ll.ll is thus not new, but the existing proofs are fairly difficult, in particular 
that of [18] which deals with the case A ^ 0. 

In the present paper we will provide alternative proofs of the energy lower 
bound and the convergence of states using the quantum de Finetti theorem 
in the same spirit as in dsiiii]. Our proofs are conceptually and technically 
simpler than those provided in [18]. The overall strategy will be explained 
in the next section. 

Our result covers the case of a rotating gas where the minimizers of the 
GP functional can develop quantized vortices. This corresponds to taking 
A{x) = HAx with 0 being the angular velocity vector. In this case, V should 
be interpreted as the trapping potential minus ^(H A x)^. The assumption 
P(a:) —>■ 00 as |x| —>• 00 is to ensure that all particles are confined to the 
system. Here our conditions on A and V are slightly more general than 
those of [12] where A is assumed to grow at most polynomially and V is 
assumed to grow at least logarithmically. 

The finite range assumption on w is not a serious restriction because we 
can always restrict the support of u; to a finite ball without changing the 
scattering length significantly. In fact, it is sufficient to assume that w is 
integrable at infinity, in which case the scattering length is well-defined. 
We can also work with a more general interaction wn ^ 0 (with scattering 
length ajv) rather than the specific choice (|1.4I) . as long as its range goes 
to zero and limisf^QQ Najsf exists; then the result in Theorem 11.11 still holds 
with a replaced by limAr_j.oo Non- In particular, if re at is chosen as in (11.61) 
for some 0 < /3 < 1, then Non -a (87r)~^ f w. The critical case /3 = 1 
considered in this paper is much more interesting because in the limit the 
true scattering length appears instead of its first order Born approximation 
(87r)-^fw. 


Acknowledgment: We thank Mathieu Lewin for helpful discussions. PTN 
is supported partially by the European Union’s Seventh Framework Pro¬ 
gramme under REA grant agreement no. 291734. NR acknowledges finan¬ 
cial support from the ANR (Mathostaq project, ANR-13-JS01-0005-01). 
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2. Overall strategy 


In this section we give an outline of the proofs of our main results, in 
order to better emphasize the key new points for the energy lower bound 
and the convergence of states. 

We shall use the following notation: Let 0 : —)■ R be a radial smooth 

Heaviside-like function, i.e. 


0 ^ 0 < 1, 6{x) = 0 for |x| ^ 1 and 6{x) = 1 for \x\ ^ 2. 

Let U : R^ —>• R be a radial smooth function supported on the annulus 
1/2 ^ |x| ^ 1 such that 


U{x) ^0 and 



U = Aira. 


For every R > 0 define 





The smooth cut-off function 9fi will be used to perform cut-offs in both 
space and momentum variables, the latter being always denoted by 


p = —iV. 

The potential Ur will be used to replace the original one. The important 
points will be that the integral of Ur yields the correct physical scattering 
length, and that we will have some freedom in choosing the range R of Ur. 

Step 1 (Dyson’s lemma). The main difficulty in dealing with the GP 
limit is that an ansatz does not give the correct energy asymptotics. In 
this regime, correlations between particles do matter, and one should rather 
think of an ansatz of the form 

N 

JJu(xi) f{xi-Xi), (2.1) 

i=l 

or a close variant, where / is linked to the two-body scattering process. We 
shall follow the approach of [18], relying on a generalization of an idea due 
to Dyson [Tj. The following lemma, proved in |19] . allows to bound our 
Hamiltonian from below by an effective one which is much less singular, but 
still encodes the scattering length of the original interaction potential. 


Lemma 2.1 (Generalized Dyson Lemma). 

For all s>0,l>e>0 and R > 2 Rq/N, we have 


Hn ^ 


N 

E 

i=i 


{hj - (1 - e)Pj6,{pj)) 


N e 


where 


N 

Wn := ^ UR{xi - Xj) 92R{xj - Xk). 


Here and in the sequel, C stands for a generic positive constant. 


( 2 . 2 ) 


(2.3) 
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Proof. Recall that the scattering length of wn is a/N. Therefore, from 
na Eq. (50) and the first estimate in (52), with {v,a,x,s) replaced by 
{wN,a/N,9s, s~^), respectively] one has 


^ N-l ^ N-l 

p^Gsip) + 2YI ~ ^ ~jr ^ ~ 

j=i j=i 


CaR^s^ 

e 


on L^(M^), for all given points yj satisfying vahij^k \yj —yk\ ^ 2i?. Since the 
left side is non-negative, we can relax the condition \yj — yk\ ^ 2i? 

by multiplying the right side with Y\k^j02R{yj — yk)- Thus for every i = 
l,2,...,iV, 


1 ^ 

pl^sipi) + 2 '^ WN{Xi - Xj) 






1 — e 

~jr 


^ Unixi - Xj) 92R{xj - Xk) - 


CaR^s^ 


j¥=i 

Multiplying both sides with 1 — e and summing over i we obtain (|2.2p . □ 


Clarification. The reader should keep in mind that we will choose R = 
R{N) —>■ 0 (actually ^ Ri$> then s ^ 00 and e ^ 0. □ 


The main point of Dyson’s lemma is that we can replace the hard inter¬ 
action potential wn by a softer one Ur which encodes the scattering length 
conveniently as f Ur = dvra. The price we have to pay for this advantage 
is twofold, however: first, we have to use all the high-momentum part of 
the kinetic energy (note that 9s{p) = 1 when p 2s); and second, the new 
potential UR{xi — Xj) comes with the cut-off “ ^k)- Together 

they really describe a “nearest neighbor” potential instead of an ordinary 
two-body potential. While the first problem is not too annoying as the low 
part of the momentum is sufficient to recover the full energy in the limit, 
the second problem is much more serious. 

Step 2 (Second moment estimate). The lower bound (12.21) leads us to 
consider the effective Hamiltonian 

HN-.= Y.hj+^-^WN (2.4) 

i=i 

where 

h := h - {1 - e)p^9s{p) - «£,*, Re,s '■= \nia[h - (1 - e)p‘^9s{p) - . 

(2.5) 


Here we use the freedom to add and remove the constant N to the 
Hamiltonian to reduce to the case h 1. In order to ensure that is 
finite, we need the extra condition 


lim 

|ir|—>-oo 


V{x) 


= 0 , 


which can be removed at a later stage, as we shall explain below. 


( 2 . 6 ) 
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We will now seek a lower bound to the ground state energy of (j2.4p . The 
philosophy, as in the previous work [18], is that if Ttv is the ground state of 
the original Hamiltonian, then roughly 

Tat « Tat f{xi - Xi) 


where / encodes the two-body scattering process and 'kAf is a ground state 
for (|2.4h . Thus the Dyson lemma allows to extract the short-range correla¬ 
tion structure and we now want to justify that 'k at can be approximated by a 
tensor power , that is, we want to justify the mean-field approximation 
at the level of the ground state of (|2.4p . 

There are two key difficulties left: 

• The effective Hamiltonian is genuinely many-body. It can be bound¬ 
ed below by a three-body Hamiltonian, but obviously one will ulti¬ 
mately have to show that the three-body contribution can be ne¬ 
glected. 

• To recover the correct energy in the limit we need to take R <C 

in order to be able to neglect the three-body contribution in the 
effective Hamiltonian. We thus still have to deal with the mean-field 
approximation in the “rare but strong collisions” limit. In other 
words, even though the effective Hamiltonian is much less singular 
than the original one, we do not have the freedom to reduce the 
singularity as much as we would like. 

It is in treating these two difficulties that our new method signihcantly 
departs from the previous works mES]. We shall rely on a strong a priori 
estimate for ground states of (|2.4p . In Lemma 13.II we assume (|2.6p and 
show that (provided R S> which is sufficient for our purpose) 

1 ^ ~ 2 

• (2.7) 

Note that a bound of this kind is not available for the original Hj\f due to 
the singularity of its interaction potential. In particular, (12.71) implies that 
every ground state Tat of satishes the strong a-priori estimate 

{'^N,hih2'^N) ^ Ce,s- ( 2 . 8 ) 

This second moment estimate is the key point in our analysis in the next 
steps. It is reminiscent of similar estimates used in the literature for the 
time-dependent problem lEiEiEniin]. 

Notation. We always denote by Ce (or Ce^g) a (generic) constant indepen¬ 
dent of s, N and R (or independent of N and R, respectively). 

Step 3 (Three-body estimate). Next we have to remove the cut-off 

6{xj - Xk) 
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in Wpi to obtain a lower bound in terms of a two-body Hamiltonian. Using 
the elementary inequality (see [THl Eq. (22)]) 

- Xfc) ^ 1 - ^ (1 - 6»2 r(xj - Xk)) 

k:k^i,j 

we have 

N 

Wn ^ ^ UR{xi - Xj) - ^ Unixi - Xj){l - 92R{xj - Xk)) (2.9) 

i^j k^i^j^k 

and we thus have only a three-body term to estimate. Since the summand 
in this term is zero except when \xi — Xj\ ^ R and \xj — Xfc| ^ 4i?, the 
last sum of (|2.9h can be removed if the probability of having three or more 
particles in a region of diameter 0{R) is small enough. This should be the 
case if R is much smaller than the average distance of particles, but 

it is rather difficult to confirm this intuition rigorously. 

In [IB], a three-body estimate was established using a subtle argument 
based on path integrals (the Trotter product formula). In this paper, we 
will follow a different, simpler approach. Instead of working directly with 
a ground state oi Hri as in [T8|, we will consider a ground state of the 
effective Hamiltonian Thanks to the second moment estimate (EID we 
can show that (see Lemma l3.4p 

N 

Ur{xi — X2)92r{x2 — Xk)'^N) ^ Ce^sNR^ (2-10) 

k=3 

The right side of (|2.10h is small with our choice ^ R. 

Step 4 (Mean-field approximation). With the cut-off in Wn removed, 
Hj\f turns into the two-body Hamiltonian 

^ _ (1 _ £)2 

Kn '■= hj -| — UR{xi — Xj) 

i=i 

for which we can validate the mean-field approximation. This is the sim¬ 
plest approximation for Bose gases where one restricts the many-body wave 
functions to the pure tensor products Since Ur converges to the delta- 

interaction with mass f Ur = dvro, we formally obtain the following approx¬ 
imation for the ground state energy 

eNL(e)'S):= inf ( (u,/lu)-|-(1 — e)^47ra / \u\'^ 

II«IIl2=i V J 

In Section 14.11 we will show that 

ini a{KN) , . 

J™ - AT -= eNL(e,s). (2.11) 

N^OO T\ 

A similar result was proved in m using a coherent state method, which 
is a generalization of the c-number substitution in [22|. In the present paper, 
we will provide an alternative proof of (12.111) using the quantum de Finetti 
theorem of Stprmer m and Hudson and Moody m- We note that this 
theorem has proved useful also in the derivation of the GP equation in 
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the dynamical case, see [1]. The following formulation is taken from m 
Corollary 2.4] (see [25] for a general discussion and more references): 


Theorem 2.2 (Quantum de Finetti). 

Let R be an arbitrary separable Hilbert space and let G with 

II^Arjl = 1. Assume that the sequence of one-particle density matrices 7 ^^^ 
converges strongly in trace class when N ^ 00 . Then, up to a subsequence, 
there exists a (unique) Borel probability measure fj, on the unit sphere SA, 
invariant under the group action of , such that 


lim Tr 

N^OO 





Vfc G N. 


( 2 . 12 ) 


This theorem validates the mean-field approximation for a large class of 
trapped Bose gases, in particular (see m and references therein) when the 
strength of the interaction is proportional to the inverse of the particle num¬ 
ber, case /3 = 0 in m- However, when the interaction becomes stronger, 
the mean-field approximation is harder to justify. The convergence (|2.11D 
with R ^ A^“2/15 pj-Qygfj [T3] by using a quantitative version of The¬ 
orem [ 22 ] valid for finite dimensional spaces [gEitts]. However, this range of 
R is too small for our purpose because we are forced to choose R <C 
in the previous steps. 

In this paper, thanks to the strong a-priori estimate (12.81) we are able to 
prove ( 12 .lip for the larger range R » As in [TTKTS] we localize the 

problem onto energy levels of the one-body Hamiltonian h lying below a 
chosen cut-off A. At fixed A, it turns out that the projected Hamiltonian is 
bounded proportionally to N. We are thus in a usual mean-field scaling if 
we are allowed to take Ai —>• 00 first, and then A —)• 00 later. Taking limits 
in this order demands a very strong control on the localization error made 
by projecting the Hamiltonian, however. This control is provided again by 
the moment estimate ( 12 . 8 p . 


Combining the arguments in steps 1-4, we can pass to the limit N ^ 00 , 
then s —>■ 00 and e ^ 0 to obtain the energy convergence CJl) under the 
extra condition (12. 6 D . In Section 14.21 we remove this technical assumption 
using a concavity argument from [18] and a binding inequality which goes 
back to an idea in |16j . 


Step 5 (Convergence of ground states). In Section [4.31 we prove the 
convergence of (approximate) ground states using the convergence of the 
ground state energy of a perturbed Hamiltonian and the Hellmann-Feynman 
principle. A similar approach was used in [18] to prove the convergence of 
the 1-particle density matrix. However, the quantum de Finetti theorem 
helps us to avoid the complicated convex analysis in [18], simplifying the 
proof significantly and giving access to higher order density matrices. 
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3. Second moment estimate 

In this section we consider the effective Hamiltonian obtained after ap¬ 
plying the generalized Dyson lemma to the original problem, namely 

^ ~ ("1 - 

Hn = Y^ hj + — Wn , 

i=i 

where h and Wn are defined in (j2.5p and (12.3p . respectively. We will work 
under the extra assumption (12.6p . Since A S M^) and V grows faster 

than at infinity, for every e > 0 we have 

{V/2 - 2e“VP)- e 

and hence 

(e/4y + y/2 - 2e-VP ^ -C*.- 

In combination with the Cauchy-Schwarz inequality, we get 

h-{l- e)p'^9s{p) ^ - 2e“VP + ^ ^ - C'e- 

Therefore, inf a{h) — 1 ^ ^ —C^ and 

^^^.(-A + H + l). (3.1) 

The key estimate in this section is the following 

Lemma 3.1 (Second moment estimate). 

Assume that dMl) holds. For every 1 > e > 0 and s > 0, if 

R = R{N) > iV^/S 

when N —>■ oo, then for N large enough we have the operator bound 

1 ^ ~ 2 

{Hn^ > . (3.2) 

i=i 

We will show in Subsection 13.31 that a convenient lower bound to Dyson’s 
potential Wn in terms of truly two-body operators follows from Lemma [3. II 
Before proving Lemma 13.11 in Subsection 13.21 we first collect some useful 
inequalities on a generic translation-invariant interaction operator W{x — y) 
that will be used throughout the paper. 

3.1. Operator inequalities for interaction potentials. We state sev¬ 
eral useful inequalities in the following lemma. In fact (13.3p is well-known 
and (|3.4I) with (5 = 0 was proved earlier in [111 Lemma 5.3]. In the sequel 
we will crucially rely on the improvement to <5 > 0, and on (13.5p which seem 
to be new. 

Lemma 3.2 (Inequalities for a repulsive interaction potential). 

For every 0 ^ IT G n L^(M^), the multiplication operator W{x — y) on 
L^((M^)^) satisfies 

0^W{x-y)^ C||1T||^3/2(k3)(-A,) (3.3) 

and, for any 0 ^ (5 < 1/4 

0i^W{x-y)^ C5||lT||ii(M3)(l - A,)i-^(l - Ay)^-\ (3.4) 
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Moreover, for all 1 > e > Q, s > Q, A ^ and 0 ^ F G 

{x-y) + W{x- y)hx 

^ -C {\\W\\l2 + (1 + 52)11^^11^3/2) (1 - A,)(l - A,). (3.5) 
Proof of Lemma \S.A We start with the 

Proof of (j3.3D . From Holder’s and Sobolev’s inequalities, we have 
(/, W{x- y)f) = jj W{x-y)\f{x, y)\‘^dxdy 

^ I {.I (/ \f(^^y)\^dx'^ 

^ C'||VF||^3/2(]k3) j (^j \Vxfix,y)\‘^dx'^ dy 

for every function / G Therefore, (13.3p follows immediately. 

Proof of (j3.4p . The estimate (13.4p with (5 = 0 was first derived in m- The 
following proof is adapted from the proof (again for (5 = 0) in |18] . Note that 
for every operator K we have K*K ^ 1 if and only if KK* ^ 1. Therefore, 
(j3.4p is equivalent to 

^W{x-y){l - - Ay)^-^y'W{x-y) ^ Cs\\Whi. (3.6) 

Let G be the Green function of (1 — A)*^”^, whose Fourier transform is given 
by 

For every function / G T^((M^)^) one has 

(/, ^W{x-y){l - A^)^-\l - Ay)^-Ww{x-y)f) 

= j f{x, y)VW(x - y)G{x - x')G{y - y’)^W{x' - y')f{x', y')dxdydx'dy’ 
^ J W{x - y)\G{x - x')\^\f{x',y')\^ + W{x' - y')G{y - y')\^\f{x,y)\'^ 

= Cs\\W\\L.{fJ) 
where 

:= I |G|2 = I |G|2 = 

which is finite for all 0 ^ (5 < 1/4. Thus ()3.6p . and hence ()3.4K . holds true. 
Simpler version of (13.5p . We are going to deduce (|3.5p from the inequality 

{-A^)W{x -y) + W{x - y){-A^) 

^ -g(||W||^3/2 + IIIFIU 2 ) (1 - A,)(l - A,). (3.7) 

By an approximation argument, one can assume that W is smooth. For 
every / G iL^(M^ X M^), a straightforward calculation using integration by 
part, and the identity S/xiW{x — y)) = —S/y{W{x — y)) gives us 

(/, ((-A,)W(x -y) + W{x - y){-Ax))f) 
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= 23^?yy V^f{x,y)\/x{W{x-y)f{x,y))dxdy 
= 2 11 \S7a,f{x,y)\‘^W{x-y) + 2‘Si jj V ^f{x,y)V ^{W{x - y))f{x,y)dxdy 
^- 23 ? JJ Vxf{x,y)Vy{W{x-y))f{x,y)dxdy 
= 2'iR ff S7y(^{Vxf{x,y))f{x,yfjW{x - y)dxdy 


= 2 ^ 


^xfix,y)Vyf{x,y) + VyVxf{x,y))f{x,y) W{x - y)dxdy. 


Using Cauchy-Schwarz and Sobolev’s inequality (13.31) . we get 




y)Vyf{x, y)W{x - y)dxdy 
|Vx/(x,y)p + \Vyf{x,y)\^ 


\W{x — y)\dxdy 


Moreover, by the Cauchy-Schwarz inequality again and (13.41) (with 5 = 0 
and W replaced by W"^), 


{VyVa;f{x,y))f{x,y)W{x - y)dxdy 


< 


\ 1/2 


\^y^xf{x,y)\^dxdy\ [jj \f{x,y)\^\W{x-y)\^dxdy\ 


\ 1/2 


^C||1U||^2(/,(1-A,)(1-A,)/). 

Thus we obtain 

(/, ((-A,)lT(x -y) + W{x - y)(-A,))/) 

^ - Ci(||W|Il3/2 + \\W\\l 2) if, (1 - A,)(l - Ay)f) 


for all f £ 


X 


. This proves (|3.7p . 


Proof of (13.5p . From the commutator relation 

PxW{x -y) = W{x- y)px + (-iVxW)(x - y) 

we find that 

{pxA{x) + A{x)px + \A{x)\'^)W{x - y) 

+ W{x- y){pxA{x) + A{x)px + \A{x)\^) 

= 2(^PxW{x - y)A{x) + A{x)W[x - y)px) + \ A{x)\^W{x - y)^ 
= 2{px + A{x))W{x - y){px + A{x)) - 2pxW{x - y)px. 

Using 

[px + A{x))W[x - y){px + A{x)) ^ 0 
and estimating pxW{x — y)px by Sobolev’s inequality (13.31) . we get 
{pxA{x) + A{x)px + \A{x)\^)W{x - y) 
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+ W{x- y){pxA{x) + A{x)px + \A{x)\‘^) 
^-C||W^||^3/2(-A,)(-A,). (3.8) 

Finally, by (I3.3p again and the Cauchy-Schwarz inequality for operators 

±{XY + Y*X*) ^ 6XX* + 6-^Y*Y, V(5 > 0, (3.9) 

we obtain 

plil - 9s{px))W{x -y) + W{x- y)plil - OsiPx)) 

> - Os{px))W{x - y)plil - OsiPx)) + 6-^W{x - y) 

^ -C\\W\\l./ 2 (5pi{l - OsiPx)? + 6-^) (-AJ 
for all h > 0. Using 1 — 9s{p) ^ l(|p| ^ 2s) and choosing 6 ~ s~‘^ gives 
Px(l - 9s{px))W{x -y) + W{x- y)plil - 9s{px)) 

^-Cs2||fU||i3/2(-A,). (3.10) 

From ()3.7I) . ()3.8p and (|3.10H . the bound (13.5p follows. □ 


3.2. Proof of Lemma 13.11 Before completing the proof of Lemma l3.ll we 
make a remark on the simpler case with the Dyson potential Wn replaced 
by a truly two-body interaction. 

Remark 3.3 (Second moment estimate with two-body interactions). 
Consider the model case 

^ _ (1 _ £')2 

Kn ■= hj H — Ufi{xi — Xj). 

1 = 1 *7^1 

By expanding and using the fact that hi ^ 0 commutes with Uji{xj — 
Xk) ^ 0 when i ^ j and i ^ k, and then using (13.51) to estimate terms of the 
form 

hiUnixi - Xj) + U{xi- Xj)hi, 

we obtain 

(3.11) 

provided that R = R{N) S> . A similar estimate also holds when h is 

replaced by the original kinetic operator h. 

We stress once again that we do not expect (|3.11l) to hold for our original 
Hamiltonian which is in the more singular regime R N-\ We thus 
need to work with the Dyson Hamiltonian, and its rather intricate nature 
makes the actual proof of Lemma 13.11 more difficult than the one we have 
sketched for (13.lip . We now proceed with this proof. 


Proof of Lemma \3.1\ We have 

N n -^2 ^ 


(R^)2 - = (l^^(W^ + WNh^) + (3.12) 

1=1 i=i 


N 
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Similarly as in Remark 13.31 the goal is to bound hiW m + W^hi from below. 
We first decompose the interaction operator as 

WN = Wa + Wb 


where 

kha = ^ UniXi-Xj) 02R{Xj-Xk), 

Wb= UR{Xi - Xj) 02R{Xj - Xk). 

i,j^2 

Estimate of Wa- By the Cauchy-Schwarz inequality (j3.9D we get 

±(hiWa + Wah) ^ N-%Wah + NWa- (3.13) 

Let us show that 

VU < (3.14) 

Indeed, for every given {xi,X 2 , ...,xn) £ the product 

Ur{xi - Xj) 02R{Xj - Xk) 
k¥^i,j 

is bounded by ||17 r||loo ^ CR~^ and it is zero except in the case 

|xi — Xj\ < R < 2R < min \xj — Xk\- 
k^^,j 

By the triangle inequality, the latter condition implies that 

Ixi — XA < R < min Ixi — xA 

k^i,j' 

and it is satisfied by at most one index j 7 ^ 1. Therefore, 

C* 

'^Ur{xi- Xj) 02R{Xj - Xk) ^ 
k^l,j 

Similarly, we have 

(j 

UR{Xi - aq) 02r{xi ~ ^k) ^ ^ 

2^2 k^l,i 

and hence (|3.14l) holds true. From (I3.13P and (I3.14p we obtain 
±(hlWa + Wahl) ^ n - Xk). 

k^i,j 

(3.15) 

Here we do not need to estimate the second term on the right side of (j3.15l) 
because this term is part of VFAr which will be controlled by in 

Estimate of Wb. We need a further decomposition 

Wb=Y. n ^2R(®i - Xk) = Wc-Wa 

i,j^2 k^i,j 
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where 


Wc := ^ iyR(Xi - Xj) 02R{Xj - Xk) 


Wd ■■= ^ UR{Xi - Xj)[l - 92R{Xj - Xl)j 02R{Xj - Xk). 

Note that 

Wc ^ 0, Wd ^ 0 and hiWc = Wchi ^ 0. 

On the other hand, by the Canchy-Schwarz inequality (|3.9I) again, 

±(hiWd + Wdhi) ^ 6hiWdhi + 5-^Wd. (3.16) 

We have two different ways to bound Wd- First, by ()3.3p and ()3.1I) . 

^1 — 92R{Xj — Xi)j < C\\l — 02i?||i3/2(l — Ai) ^ CeR^hi. 

Since here i,j ^ 2, both sides of the latter estimate commute with 
UR{Xi-Xj) 92R{Xj-Xk), 

and we deduce that 

(l - 02R{Xj - Xl)^ UR{Xi - Xj) 62 R{Xj - Xk) 

^ CeR^hiUR{Xi - Xj) 92R{Xj - Xk). 

Taking the sum over i,j ^ 2 we obtain 

Wd ^ CcR%Wc. 

Second, let us show that 

C 




R^' 


(3.17) 

(3.18) 


Indeed, for every given (xi,X 2 , ...,xn) £ (K^)-^, the product 

UR{Xi - Xj) (^1 - 02R{Xj - Xl)^ 92R{Xj - Xk) 


ky^i,i,j 


is zero except in the case 


|xj — Xj\ < R, \xj — xi| < 4ii, min \xj — Xk\ > 2R. (3.19) 

By the triangle inequality, (13.191) implies that the ball B{xi,5R) contains 
B{xi, R/2), B{xj, R/2), and the balls B{xi, R/2), B{xj,R/2) do not inter¬ 
sect with B{xk, R/2) for all k ^ Since B{xi,5R) can contain only a 

finite number of disjoint balls of radius R/2, we see that there are only a 
finite number of pairs {i,j) satisfying (13.191) . Thus we can conclude that 

Wd ^ C\\Ur\\l<^ ^ CR-\ 

From ()3.16l) . ()3.17l) and (I3.18p . we obtain 

hiWi + Wjii = hiWd + Wdhi + 2hiWc 
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Choosing 6 I'N-' we get 

hiWb + Wbhi ^ (3.20) 

rC 

Conclusion. From (j3.15p and (j3.20p . we get 

hi Wat + WnHi ^ (hi)^ 

-2iV ^ Unixi-Xj) 02Rixj-Xk). 

Summing the similar estimates with 1 replaced by i and using 
N 

Unixi - Xj) 62 R{xj - Xk) = 2 Wn 

^=1 kj^i,j 

we find that 

/ C C \ ^ ~ 

[heWN + WNhe) ^ - i ^ ^ - 2NWn. 

i=i ^ ^ 1=1 


Therefore, coming back to (I3.12p we conclude that (completing a square in 
the last inequality) 


N 


("jv)'- E* 


i=i 


(l-£) 

N 


E 

1=1 


+ 


+ iV2 


- 


- 


(]^|s + ll) E(h)^ - 2(1 - 

^ ^ 1=1 

+ E(h)^ - IV^- 


N'^R^ NR 


i=i 


When R^ N we have 


C Cr 


+ < 1 


N‘^R^ NR 


and hence 


N 


{Huf ^ 2 h,hj + (1 - o(l)) J2{hif - N^ 


£=1 


which yields the result, recalling that in our convention h ^ 1. □ 


3.3. Three-body estimate. A first consequence of the second moment 
estimate in Lemma 13.11 is that we can conveniently bound Dyson’s Hamil¬ 
tonian from below by a two-body Hamiltonian. This is done by first using 
a simple bound in terms of a three-body Hamiltonian, and then bounding 
the unwanted three-body part. 
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Lemma 3.4 (Three-body estimate). 

Assume the extra condition (USD holds. For every 1 > e > 0 and s > 0, if 
R = R{N) > A^-2/3, then 

^ Unixi - Xj) ^ (1 - 02 R{xj - Xk)) ^ Ce,s^{HN)^. (3.21) 

Consequently, 

Hn>Y1 + Sv ^ ~ ~ (3-22) 

1 = 1 * 7^1 

Note the error term involving (/^tv)^) which is well under control since we 
are interested in its expectation value in a ground state. 

Proof. By (|3.3p and ()3.ip we have 

{'i- - d 2 R{x 2 - Xk)) ^ Ce^sF?lvk for fc ^ 3. 

Since Ur{xi — X 2 ) commutes with both sides, we get 

Ur{xi - X2) ^(1 - 62r{x 2 - Xk)) ^ Ce,sF?UR{xi - X 2 ) ^^3 

= (Hn -hi-h2-{l- efN-^WN) Ur{xi - X2) 

+ ]^Ce,sR^UR{xi - X2) (Hm -hi-h2-{l- e)^N-^WN) 

^ ]^Ce,sR^ {HnUr{xi - X2) + Ur{xi - X2)Hn^ 

+ -Ce,sR^ ^ {hjUR{xi - X2) + Ur{xi - X2)hj^ . (3.23) 

1=1 

In the last estimate we have used Wn ^ 0. Thanks to (|3.5I) and (121]) , we 
get for all j = 1, 2, 

hjUR{xi - X 2 ) + Ur{xi - X2)hj 

^ -C,,sR-^/^{l - Ai)(l - A 2 ) ^ -Cs,sR-^^%h2. (3.24) 

On the other hand, by the Cauchy-Schwarz inequality (13.9p and (|3.4I) (with 
(5 = 0 and IT = Ur) 

HnUr{xi - X2) + Ur{xi - X2)Hn 
^ 6HnUr{xi - X2)Hn + 6~^Ur{xi - X2) 

^ Cs,sdHNhih2Hi\f + C'e,s5“^/ii/i2 (3.25) 

for all 5 > 0. Choosing 6 = N~^ and using R~^/‘^ ^ N, we deduce from 
p3.23p . (I3.24P and p3.25l) that 

Ur{xi - X2) ^(1 - 62r{x 2 - Xk)) ^ Ce,sR^{N~^HNhJl2HN + Nhjl2^. 
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By symmetrization with respect to the indices, we find that 
^ Ur{xi - X 2 ) ^(1 - ^2R(Xj - Xk)) 

^ Ce,sR^ (n-^Hn Y, hihjHN + nY • 

i¥=j 

Combining with the second moment estimate (|3.2H we obtain (I3.2ip . From 
the three-body estimate (I3.2ip and the elementary inequality (j2.9h . the op¬ 
erator bound (j3.22h follows. □ 

4. Energy lower bound and convergence of states 

4.1. Mean-field approximation. We are now reduced to justifying the 
mean-field approximation for a new Hamiltonian with the two-body inter¬ 
action Ur{x — y) which converges to a Dirac delta much slower than the 
original one. The analysis in this section provides an alternative to the 
coherent states method of [18] . 

Proposition 4.1 (Mean-field approximation). 

Assume that (j2.6p holds. For every 1 > e > 0 and s > 0, if 

> ii = R{N) > iV“2/3 

then 

lim -- f = inf (^(u,/itt)-|-(1 — e)^47ra [ =:eNL(£ 5 s). 

N \\u\\^2=l V J J 

(4.1) 

Proof. The upper bound in (14.ip can be obtained easily using trial states of 
the form . For the lower bound, let us consider a ground state Tw of 
Hjsf (which exists because h has compact resolvent). Using the ground state 
equation, we find that 

(Tw, = (inf aiHN))’^ ^ {Ce,sN)’^ (4.2) 

for all A; € N. In particular, the second moment estimate (13.2D implies that 

N, hih2'it n) ^ C'e.s (4.3) 

and the operator estimate (j3.22D implies that 

lim inf ^ W (1 _ £)2 . 

N—^00 JV N—^00 V V '^nJ V ^nJJ 

(4.4) 

(k) ^ 

Here 'j- ' is the A:-particle density matrices of Tw and Ur is understood as 

'Ejv 

the multiplication operator Ur{x — y) on . Since Tr ^ is bounded 

uniformly in N and h has compact resolvent, up to a subsequence we can 
assume that 7 -^ converges strongly in trace class. By the quantum de 
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Finetti Theorem 12.21 up to a subsequence we can find a Borel probability 
measure Jl on the nnit sphere SSj snch that 


lim Tr 7 -^^ 

N^oo '^N 


/' 




= 0, Vfe G N. 


We will show that 

Itoinf (tv + (1 - e)^TV (C/OT?()) 

^ /iti) + (1 — 6)^4™ / 


(4.5) 


(4.6) 


and then the lower bonnd in (Bm follows immediately. Since h is positive 
and independent of N, (14.5p and Fatou’s lemma imply 


liminf Tr ^ / {u,hu)dfj,(u). 

IV-voo V Viv/ J 


It remains to prove 


lim inf Tr 

N^OO 



(4.7) 


(4.8) 


Note that (j4.8p does not follow from ()4.5p and Fatou’s lemma easily because 
Ur depends on i? = R{N). We need to replace Ur by an operator bounded 
independently of N. Since h has compact resolvent, for every A ^ 1 the 
projection 

Pa := 1(^ ^ A) 

has finite rank. Let us denote 


n:= 

Since Ur ^ 0, we can apply the Cauchy-Schwarz inequality (18.9p with 
X = Pf^u]!"^ and Y = 

to obtain 


Ur = {Pf +Ii)UR{Pf^ Pli) 

= PfURPf + UUrU + PfURU + UURPf 
^ Pf‘^URPf‘^ - S-^UUrU - SP^^UrP^'^ 


for all d > 0. Using the operator bound 
power is operator monotone [3] we have 


and the fact that the 4/5-th 


Ur{xi - X2) ^ C'||Ur||^i(1 - Ai) 4/5(1 - ^ C,,,(/ii)"/5(/i2)"/^ 


(4.9) 


Therefore, 

P^^UrP^ ^ Ce,shji 2 and UUrU C,,sA-^/%h 2 . 

Here in the second estimate we have used — Pa ^ A“^/®(/i)^/^, which is 
a consequence of the definition of Pa- Thus 

Ur - P^^UrP^ ^ + 6A-^/^)hJi2. 
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If we choose 6 = A and take the trace against 7 -^^ , then by the a-priori 
estimate (j4.3p we find 

Tr - Tr ^ -C,,sA-^/^^. (4.10) 

On the other hand, from dMD and the definition of Pa; it follows that 
the operator norm \\P^‘^UrP^‘^\\ is bounded uniformly in N for fixed A. 
Therefore, the strong convergence (14.5p implies that 

Jim^ (^TV (Pf - I {{PAu)^\UR{PAu)^^)dJi{u)^ =0. 

(4.11) 


Since the left side of (j4.7p is finite, every function u in the support of djl 
belongs to the quadratic form domain Q{h) oih and hence Pau —>■ u strongly 
in Q{h). Using the continuous embeddings Q{h) C C L^, we get 


lim lim((PAu)^^,PR(PAu)^^} = lim WPAuWii = 
A—^(50 R —^0 A—^00 


I il4 


By Fatou’s lemma, 

liminfliminf [{(Pau)^"^,UR{PA u)^‘^)d'il{u) ^ Aira [ Wulliidjliu). (4.12) 

A^oo N^oo J J 

The desired convergence (14.8p follows from p4.inp . p4.1ip and (j4.12K . □ 


Remark 4.2 (Mean-field approximation with two-body interactions). 

From the preceding proposition we obtain easily the convergence (12.111) men¬ 
tioned in Section 0 because Hj\[ ^ Kja- In fact, Kja satisfies the second 
moment estimate \3.11\) (cf. Remark \3.3i) . and hence (12.111) can be proved 
directly. In particular, the method can be used to derive the energy asymp¬ 
totics when the interaction potential is given by 11.51) .' for [3 < 2/3, Step 1 
(and thus also Step 3) in the proof are not needed. One can also obtain some 
explicit error estimate in Proposition |/.i| and (12.lip by using a quantitative 
version of the quantum de Finetti theorem as in m Lemma 3.4]. 


4.2. Convergence of gronnd state energy. We now conclude the proof 
of the convergence of the ground state energy. There are two things left to 
do: remove the high momentum cut-off in the final effective functional, and 
relax the additional assumption ( 12 . 6 D . 


Proof of energy convergence (11.9p . The upper bound in (11.9p was proved 
in [26]. The proof of the lower bound is divided into three steps. 

Step 1. We start with the simple case when the extra condition (12.61) holds 
true. Recall that we are choosing 

> P = R{N) > A^“2/3_ 

From Lemma ED and Proposition 14.11 it follows that for every 1 > e > 0 
and s > 0 , 


N^oo N 


^ lim inf 
N^oo 


(inf a(PAr) ^ ^ 

[ N ^ J 


eNL(e, s) + He,s- 
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Thus to obtain the lower bound in (|1.9I) . it remains to show that 

lim lim (cnlCe, s) + Ke,s) = ecp- (4-13) 

s^oo 

The upper bound in (14.131) is trivial as £’NL(^f) + He,s ^ ^^gp(^^)- The lower 
bound in (j4.13l) can be done by a standard compactness argument provided 
in [18]. We recall this here for the reader’s convenience. Let be a ground 
state for eNL(e)'S), namely 


eNL(e, s) = {us,s, hue^s) + (1 - e)^47ra j \us,s 


From (13.11) it follows that (ue,s, (—A + V)u£^s) is bounded uniformly in s. 
Since —A + V has compact resolvent, for every given e > 0 there exists a 
subsequence sj —)• oo such that converges strongly in and pointwise 
(in both p-space and x-space) to a function u^. By Fatou’s lemma we have 

liminf / \us s ■{x)\‘^dx ^ / \us{x)\‘^dx, 
j^oo J ' ^ J 

liminf j p^il - Osj{p))\uE,sj{p)\‘^dp ^ J p‘^\ue{p)\‘^dp. 

Next, using (12.61) as before we have 

ep^ + pA + Ap + \ A\‘^ + V + Ce^0 

for some Cg ^ 0. Using Fatou’s lemma again and the strong convergence in 
we deduce 

Inn inf (^Ue,sj, +pA + Ap+ + U + 

^ (ue, {ep^ +pA + Ap+ + v'^Ue^ . 
Combining these estimates, we get 

Inninf (^eNL(e, Sj) + > {ue, hue) + (1 - e)^47ra j \ue\‘^ ^ (1 - e)^eGp. 

Taking e —>■ 0 we obtain the lower bound in (|4.13l) . 

Step 2. From now on we do not assume (12.6|) . Let us introduce the Hamil¬ 
tonian 

M 

Hm,n ■= '^hj+ ^ WN{xi - Xj) 

j=l l^i<j^M 

and denote E{M, N) its (bosonic) ground state energy. In this step we will 
prove the lower bound in (11.91) using the additional assumption 

E{N,N) - E{N-1,N) ^C. (4.14) 

We will find a function / : —)■ M_|_ growing faster than |T|, namely 


l.m 

hKoo J\x) 

such that for a ground state Tjv for Ff at we have 

(Tv, /^(a;i)Tv) ^ C. 


(4.15) 


(4.16) 
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Once this is achieved we get 

N ^ N ^ 

for every 77 > 0. Since the growth condition ( 12 . 6 p holds true with V replaced 
hy V + 77 /^, we can apply the result in Step 1 to the Hamiltonian 

N 

HN + r]^f{xj) 
j=i 

for every given 77 > 0. Then the lower bound in (11.911 follows by taking 
77 —)• 0 . 

Now we find such a function /. We will establish a simple binding in¬ 
equality using an idea in m- From the ground state equation = 

E{N,N)^i\f, it follows that 

E{N, N){^n, /'(xiv)Tjv) = f{xN)HN,N^N). (4.17) 

By the variational principle and (I4.14p we have 

Hn,n -hN> i77V-i,A ^ E{N - 1, iV) ^ E{N, N) - C. 

Note that f'^{xN) commutes with all terms in the latter inequality. If / is 
bounded and sufficiently regular, we have the IMS-type formula 

^ifh + hf) = fhf - |V/|2 ^ vf - |V/|2, (4.18) 

and we deduce from (j4.17p that 

(Ttv, (v{xN)f{xN) - |V/(x7v)|2 - Cf{xN))^N) ^ 0. (4.19) 

Note that if we choose f{x) = for some constant 6 > 0, then (14.1511 
follows from the assumption (11.211 . Moreover, heuristically (|4.16ll follows 
from (j4.19ll as Vf'^ grows faster than |V/p -|- C/^. To make this idea 
rigorous, let us apply (I4.19p with f{x) replaced by 

grix) = exp [b [r - ||x| - r\]^] . 

Note that gr{x) = when |a;| ^ r and gr{x) = 1 when |x| ^ 2r. We can 
thus apply (I4.18p to gr- 
Moreover, 

Vgl - \Vgr? - Cgl ^{V-b^- C)gl 

^ gl -{b^ + C + l)gll {V^b^ + C + l) 

^ 9r ~ Co 

for some constant Cq independent of r > 0. Here we have used the fact that 
5 ^ 1(17 ^ 6 ^ -|- C -|- 1 ) is bounded independently of r > 0 , which follows from 
the assumption lim| 2 ,|^cxD I7(x) = -|-oo. Thus (I4.19p gives us 

{'^N,grixN)'^N) ^ Co 

for all r > 0. Taking r —>■ oo we obtain (I4.16p with /(x) = 


Step 3. Now we explain how to remove the additional assumption (I4.14p . 
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This can be done by following the strategy in [18], which we recall quickly 
below for the reader’s convenience. 

By choosing trial states , we get the upper bound 

E{N, N) ^ CoN 

for some constant Cq > 2eQp. For every € N, we denote by M = M{N) 
the largest integer ^ N such that 

E{M{N),N) - E{M{N) - 1, iV) ^ Cq. (4.20) 

Then by the choice of M{N) we obtain 

E{N, N) - E{M{N),N) ^ {N - M{N))Co. (4.21) 


We can find a subsequence Nj —)• oo such that M{Nj)/Nj —)• A € [0,1]. 
Since (14.201) holds with M = M{Nj) we can apply the result in Step 2 with 
w replaced by Xw and hnd that 

liminf ^ eGp(Aa) ^ AeGp(a). (4.22) 

j—)-oo JVij 

Here eGp(Aa) is the Gross-Pitaevskii energy with a replaced by Xa and the 
last inequality in (|4.22l) is obtained by simply ignoring part of the one-body 
energy in the corresponding Gross-Pitaevskii functional. From (|4.2ip and 
(|4.22l) . it follows that 


eGp(a) ^ liminf 
j^oo 


EiNj,N,] 

N. 


^ liminf 
j^oo 

v2 


fE{M{N,),J^ ^ - M{N,) 




N, 


Ni 


^ A eGp(a) + Cq{1 — X). 


Since 

eGp(a) ^ A^eGp(a) + 2(1 — A)eGp(a) 

and Cq > 2eGp(o), we must have A = 1. Thus M{N)/N —)• 1 for the whole 
sequence and 


lim inf 

N^OO 


E{N, N) 
N 


= lim inf 
j-i>oo 


E{Nj,N,] 

N, 


> eGp(a). 


This completes the proof of the energy convergence (jl.Op . 


□ 


4.3. Convergence of density matrices. Now we prove the convergence 
of ground states in (|1.10l) by means of the Hellmann-Feynman principle. For 
V € L^(M^) and G N we will perturb EIn hy 

■■= 7 ^ E .ir 

Here denotes the operator acting on the £-body 

Hilbert space of the ii-th,..., i^-th variables. We have the following extension 

of (frgp . 

Lemma 4.3 (Energy lower bound for perturbed Hamiltonians). 

We assume dO]), dop and dni). For every v G L^(M^) and i we have 

liminf ^ jnf ('£’gp(m) - \ {v,u)\A . (4.23) 

N—^oo Jy ||u||^ 2 =l 2 
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Proof. We first work under the extra condition (12.61) . and then explain how 
to remove it at the end. Let 1 > e > 0 and s > 0 and 

> i? = R{N) > iV2/3_ 

Recall that from (j2.2l) we have 

Hn - >Hn- + NKe,s - Ce,sNR\ (4.24) 


Let 4>7 v be a ground state for Since is bounded uniformly 

in N, (14.2p still holds true with 'Ljv replaced by <hAr, namely 

($JV, ^ {Ce,sNf (4.25) 


for all A: € N. Combining (I4.25P with the three-body estimate in Lemma [3.4l 
we get the following analogue of (14.41) 


. inia{HN 
hm ml-— 

N—>-oo N 



= liminf 

N^OO 


(<l>Ar, {Hn - 

N 


(4.26) 


^ liminf 

TV— >-CX) 


(K'O + (1 - (^«C) “ (l■“®'>('’®L£)) 


Moreover, (|4.25p and the second moment estimate (13.21) imply the a-priori 
estimate {^N-,hih 2 ^n) ^ C'e.s- Therefore, we can estimate the right side 
of (j4.26l) by proceeding exactly as in the proof of Proposition 14.11 More 
precisely, by the quantum de Finetti Theorem 12.21 we can find a Borel 
probability measure /r$ on the unit sphere SS^ such that, up a subsequence. 


lim Tr 

TV^OO 



/' 




= 0 , 


Vfe e N. 


Using (|4.6I) with 'Ltv replaced by <hAr and employing the fact that \v®^){v^ 
is bounded, we obtain 


lim inf (Tr 

TV^OO 


( 1 ) 

'J'jv 


( 2 ) 

'J'jV 






+ (l-e)^Tr 

J (^{u,hu) + {1 — e)‘^4:7ra j \u\'^ — \{v,u)\^^'^d^q>{u). (4.27) 


From ()4.24l) . ()4.26l) and (I4.27p . it follows that 
inf a{HN - S^^i) 


lim inf ■ 

TV—>-oo 

^ inf 
Il“lli2=l 


N 


■/ 


{u,hu) + {\ — eYAira I — |(u, u)p^ ) -|- 


The lower bound p4.23l) follows by passing to the limits s —)■ 0 and then 
e —)• 0 as in the proof of (I4.13p . 

To remove the assumption (j2.6l) we may use the argument in Subsection 
1121 The only extra difficulty is that when dealing with the analogue of 
P4.17I) with Hm,n replaced by — Sy^, we have to take care of the 

operator f‘^\v){v\ = |/^u)(u| which may be unbounded as f{x) = with 
6 > 0 and v is merely in L^(]R^). However, we can still proceed with all 
functions v in L^(M^) which have compact support. Then after obtaining 
the lower bound (j4.23l) with those nice functions v, we can extend the lower 
bound to all functions v in L^(M^) by a standard density argument. □ 
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Now we are able to prove the convergence of density matrices. 

Proof of state convergence (ll.lOp . Let 'Lat be an approximate ground state 
for H]\f as in Theorem ll.il For every v G L^(R^) and £ G N, from the upper 
bound in (jl.9l) and the lower bound in Lemma 14.31 we have 


limsupTr 

N->-oo ^ ^ 


1 - f N, Hisf^ n) - Sy £)'^N) 

limsup I -—-—- 


N—>-oo \ 


N 


N 




€ 


f inf a{HN) inf a{HN - 

'7^1 [— -s-j 

ecp - inf (TGp(n) - |(n,n)|^^) . 
||«ih2=l ^ ^ 


Here v is not necessarily normalized. Therefore, we can replace v by 
with A > 0 and obtain 

limsupTr ^ j (^ecp - (^Gp(^^) - AKn,n)p^)^ . 

(4.28) 

With given v and i, for every A > 0 let ux be a (normalized) minimizer for 
u i-A £’Gp(n) — \\{v,u)\^^. Since {ux,hux) is bounded and h has compact 
resolvent, there exists a subsequence Xj —>■ 0 such that ux. converges to uq 
in L^. By Fatou’s lemma, uq is a minimizer of £’gp(^)- Moreover, 

lim sup( eGP - „ inf (fGp(w) - Aj|(n, ) 
j^OO V IfIL2-1 ^ ^ J 

^ limsup^ (fGp(MA^) - (SQniuXj) - Aj|(ii,'WA,)P^)) = \{v,uo)\^^. 

(4.29) 

From ()4.28l) and (I4.29p . we conclude that for every v G L^(M^) and £ € N, 
limsupTr ^ sup \{v,u)\^^ (4.30) 

TV—>^oo ^ ^ 11 ^ \/1/-i T~, 


uGMqp 


where AIgp is the set of minimizers of £’Gp(n). 

Note that also in [18] the upper bound (I4.30p with i = 1 was proved, and 
from it the convergence of the one-particle density matrices was deduced 
using an abstract argument of convex analysis. In the following, we will 
provide a simpler way to conclude the convergence of density matrices from 
(I4.30p . using the quantum de Finetti Theorem 12.21 Indeed, by Theorem 
12.21 as before, up to a subsequence of Ttv, there exists a Borel probability 
measure /r on the unit sphere SSj such that 


lim Tr 
N—^OO 


itl- j 


= 0, \/kG N. 


(4.31) 


We will show that ^ is supported on AIgp- From (|4.30p and ()4.3ip . we get 
\{v,u)\‘^’^dfj,{u) ^ sup |(i;,tt)p^, Vi; G L^(M^), A: G N. (4.32) 


/ 
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We assume for contradiction that there exists vq in the support of fx and 
uo ^ AIgp- We claim that we could then hnd 5 G (0,1/2) such that 

|(u,u)| ^ 1 — 3(5^, \/u e M.gp,'^v £ B (4.33) 

where B is the set of all points in the support of // within a L^-distance 
less than <5 from vq. Indeed, if that were not the case, we would have two 
sequences strongly converging in 


Vn Vo, Un^Uo£ MqP 


with \\un — Vn\\ —>■ 0, and thus vq £ M.gp- Here we have used that AIgp is 
a compact subset of L^(M^). 

On the other hand, by the triangle inequality. 


|(u,u)| > 


IttlP + llulP — \\u — ulP 


^ 1 — 26^, Vu, V £ B. 


(4.34) 


Combining (|4.32H , (|4.33l) and (|4.34l) we find that 


{n{B)f{l-25‘^f^f^ [ j \{v,u)\^^d^i{u)d^i{v) 

Jb Jb 

^ f sup \{v,u)\'^^dfi{v) ^ (4.35) 

JB mGAIgp 


for all A: G N and hence, taking A: —>■ oo, //(H) = 0. However, it is a 
contradiction to the fact that vq belongs to the support of // and // is a 
Borel measure. Thus we conclude that fi is supported on AIgp and the 
proof is complete. □ 
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